Evidence for the pair- breaking process in li61i/ Sn 
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The nuclear level densities of 116 ' 117 Sn below the neutron separation energy have been determined experi- 
mentally from the ( 3 He,o:y) and ( 3 He, 3 He'y) reactions, respectively. The level densities show a characteristic 
exponential increase and a difference in magnitude due to the odd-even effect of the nuclear systems. In addi- 
tion, the level densities display pronounced step-like structures that are interpreted as signatures of subsequent 
breaking of nucleon pairs. 

PACS numbers: 21.10.Ma, 24.10.Pa, 25.55.-e, 27.60.4-j 



I. INTRODUCTION 

Nuclear level densities are important for many aspects of 
fundamental and applied nuclear physics, including calcula- 
tions of nuclear reaction cross sections. The level density of 
excited nuclei is an average quantity that is defined as the 
number of levels per unit energy. The majority of data for 
nuclear level densities are obtained from two energy regions. 
At low excitation energy, the level density is obtained directly 
from counting of low-lying levels [1]. As the excitation en- 
ergy increases, the level density becomes large and individual 
levels are often not resolved in experiments. Therefore, the 
direct counting method becomes impossible. Nuclear reso- 
nances at or above the nucleon binding energy provide another 
source of level density data [2]. Between these two excitation 
energy regions, the level density is often interpolated using 
phenomenological formulas (j, 0, |5|]. It is in this energy re- 
gion the present measurements focus. 

Recently, an extension of the sequential extraction method, 
now referred to as the Oslo method, was developed by the 
Oslo Cyclotron Group. The Oslo method permits a simul- 
taneous determination of the level density and the radiative 
strength function 10,0]. For both of these quantities, the ex- 
perimental results cover an energy region where there is little 
information available and data are difficult to obtain. How- 
ever, the limitation of this method is that the results must be 
normalized to existing data - from the discrete levels and neu- 
tron resonance spacings for the level density and to the total 
radiative width for the radiative strength function. Thus, the 
new and main achievement of the Oslo method is to estab- 
lish the functional form of the level density and the radiative 
strength function in the above specified energy region. 

In this work we present results for the level density in 
1 16,1 i7 Sn for tne exc i tat i on energy < E < S n - 1 MeV. The 
radiative strength functions of 116117 Sn have been published 
elsewhere [8J]. The experimental set-up is briefly described in 
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Section II, followed by a discussion of the analysis and nor- 
malization procedure. The experimental results for the level 
density are given in Section III, and the determination of var- 
ious thermodynamic quantities are presented in Section IV. 
Conclusions are drawn in Section V. 



II. EXPERIMENTAL PROCEDURE AND DATA ANALYSIS 

The experiment was carried out at the Oslo Cyclotron Labo- 
ratory (OCL) using a 38-MeV 3 He beam. The self-supporting 
117 Sn target had a thickness of 1.9 mg/cm 2 . The reaction 
channels 117 Sn( 3 He,ay) 116 Sn and 117 Sn( 3 He, 3 He'y) 117 Sn 
were studied. 

The experiment ran for about 11 days with an average 
beam current of k 1.5 nA. Particle-y coincidence events 
were detected using the CACTUS multidetector array. The 
charged particles were measured with eight Si particle tele- 
scopes placed at 45° with respect to the beam direction. Each 
telescope consists of a front Si AE detector with thickness 140 
jim and a back Si(Li) E detector with thickness 3000 jJ,m. An 
array of 28 collimated Nal y-ray detectors with a solid-angle 
coverage of w 15% of An was used. In addition, one Ge detec- 
tor was used in order to estimate the spin distribution and de- 
termine the selectivity of the reaction. The typical spin range 
is / - 2 - 6H. 

Figure Q] shows the singles a-particle spectrum (upper 
panel) and the a-y coincidence spectrum (lower panel) for 
116 Sn. The two peaks denoted by + and 2 + are the transfer 
peaks to the ground state and the first excited state, respec- 
tively. The strong transfer peak at E = 3.2 MeV is composed 
of many states found to be the result of pick-up of high- j neu- 
trons from the g 7 / 2 and h n/ / 2 orbitals 0. Another strong 
transfer peak centered around E = 8.0 MeV is new and may 
indicate the neutron pick-up from the gg/2 orbital. The counts 
in the coincidence spectrum decrease for excitation energies 
higher than S„ due to lower y-ray multiplicity when the neigh- 
boring A — 1 isotope is populated at low excitation energy. 

The particle-y-ray coincidence spectra were unfolded using 
the response functions for the CACTUS detector array fiUl . 
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FIG. 1: a) Singles a-particle spectrum, and b) a-J coincidence spec- 
trum from the 1 17 Sn( 3 He, a) 1 16 Sn reaction. 



The first-generation matrix (the primary-y-ray matrix) con- 
tains only the first y-rays emitted from a given excitation en- 
ergy bin; this matrix is obtained by a subtraction procedure de- 
scribed in Ref. B12I1 . This procedure is justified if the y decay 
from any excitation energy bin is independent of the method 
of formation - either directly by the nuclear reaction or indi- 
rectly by y decay from higher lying states following the initial 
reaction. This assumption is clearly valid if the same states 
are populated via the two processes, since y-decay branching 
ratios are properties of levels. When different states are popu- 
lated, the assumption may not hold. However, much evidence 
suggests ]7|] that statistical y decay is governed only by the 
y-ray energy and the density of states at the final energy. 

In the data analysis, the particle- y-ray coincidence matrix is 
prepared and the particle energy is transformed into excitation 
energy using the reaction kinematics. The rows of the coinci- 
dence matrix correspond to the excitation energy in the resid- 
ual nucleus, while the columns correspond to the y-ray energy. 
All the y-ray spectra for various initial excitation energies are 
unfolded using the known response functions of the CACTUS 
array II in . The first-generation y-ray spectra, which consist 
of primary y-ray transitions, are then obtained for each exci- 
tation energy bin lfl2ll . 

The first-generation y-ray spectra for all excitation energies 
form a matrix P, hereafter referred to as the first-generation 



matrix. The entries of the first-generation matrix P are the 
probabilities P(E,Ey) that a y-ray of energy Ey is emitted 
from an excitation energy E. This matrix is the basis for 
the simultaneous extraction of the radiative strength function 
and the level density. According to the Brink- Axel hypothe- 
sis lfl3l [Till , a giant electric dipole resonance can be built on 
every excited state with the same properties as the one built on 
the ground state, that is, the radiative strength function is inde- 
pendent of excitation energy and thus of temperature. Many 
theoretical models do include a temperature dependence of the 
radiative strength function lfl5l[l6ll . However, the temperature 
dependence is weak and the temperature change in the energy 
region under consideration here is rather small. Therefore, the 
temperature dependence is neglected in the Oslo method. 

The first-generation matrix is factorized into the radiative 
transmission coefficient, which is dependent only on Ey, and 
on the level density, which is a function of the excitation en- 
ergy of the final states E — Ey\ 



P(E,Ey)°< £r(E r )p(E-E r ). 



(1) 



The functions p and & are obtained iteratively by a globalized 
fitting procedure [6]. The goal of the iteration is to determine 
these two functions at ^ TV energy values each; the product of 
the two functions is known at ~ N 2 /2 data points contained 
in the first-generation matrix. The globalized fitting to the 
data points determines the functional form for p and ST . The 
results must be normalized because the entries of the matrix P 
are invariant under the transformation [6] 

p{E-E r ) = Aexp[a(E-Ey)]p{E-Ey), (2) 
#(E r ) = Bexp(aEy)£?(Ey). (3) 

In the final step, the transformation parameters A, B, and a 
which correspond to the most physical solution must be deter- 
mined. Details of the normalization procedure are described 
in several papers reporting the results of the Oslo method, see 
for example lfl7i [Tol . In the following, we will focus on the 
level density and thermodynamic properties. 



III. LEVEL DENSITIES 

The coefficients A and a relevant for the nuclear level den- 
sity are determined from normalizing the level density to the 
low-lying discrete levels and the neutron resonance spacings 
just above the neutron separation energy S„. For 117 Sn, we 
used s- and p-wave resonance level spacings (Do, D\) taken 
from |01 to calculate the total level density at S„ (see Ref. [6] 
for more details on the calculation). Since there is no exper- 
imental information about Do or D\ for 116 Sn, we estimated 
p(S n ) for this nucleus based on systematics for the other tin 
isotopes |2j,|4(], and assuming an uncertainty of 50%. The level 
spacings and the final values for p(S„) are given in Table U 

The experimental level density p is determined from the 
nuclear ground state up to <~ S„ — 1 MeV. Therefore, an in- 
terpolation is required between the present experimental data 
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TABLE I: Parameters used for the back-shifted Fermi gas level density and the calculation of p(S n ). 



Nucleus £p a i r C\ a Dq 




a(S„) S n 


P(S„) T) 


(MeV) (MeV) (MeV^ 1 ) (eV) 


(eV) 


(MeV) (10 5 MeV" 


') 


llb Sn 2.25 -1.44 13.13 




4.96 9.563 


4.12(206) 0.46 


I17 Sn 0.99 -1.44 13.23 507(60) 


155(6) 


4.44 6.944 


0.86(25) 0.40 



and p evaluated at S„. The back-shifted Fermi gas level den- 
sity with the global parameterization of von Egidy et aJ. H, 



p{E) = n 



exp(2-\ 



(4) <T 



is employed for the interpolation 1 . The intrinsic excitation en- 
ergy is given by U = E — £p a j r — Ci, where C\ = — 6.6A~ 0,32 
MeV is the back-shift parameter and A is the mass number. 
The pairing energy £paix is based on pairing gap parameters 
A p and A n evaluated from even-odd mass differences 11911 ac- 
cording to [20]. The level-density parameter a and the spin- 
cutoff parameter a are given by a — 0.2 1A 087 MeV -1 and 
a 1 = 0.0888A 2//3 a7\ r espec tively. The nuclear temperature T 
is described by T = y/V /a MeV. The constant 77 is a parame- 
ter applied to ensure that the Fermi gas level density coincides 
with the neutron resonance data. All parameters employed for 
the Fermi gas level density are listed in Table [I] 

Figure |2] shows the normalized level densities of 116 Sn and 
117 Sn. The full squares represent the results from the present 
work. The data points between the arrows are used for nor- 
malizing to the level density obtained from counting discrete 
levels (solid line) and the level density calculated from the 
neutron resonance spacing (open square). The discrete level 
scheme is seen to be complete up to excitation energy » 3.5 
MeV in 116 Sn beyond which the level density obtained from 
discrete levels starts to drop. For 117 Sn the discrete level den- 
sity is complete only up to w 1.5 MeV. The new data of this 
work thus fill the gap between the discrete region and the cal- 
culated level density at S n . 

From Fig. |2j we observe that pronounced step-like struc- 
tures are present in the level densities of both 116117 Sn. In the 
following section, nuclear thermodynamic properties are ex- 
tracted using the present level density results, and these struc- 
tures are investigated in detail. 



IV. THERMODYNAMIC PROPERTIES 

The entropy S(E) is a measure of the number of ways to ar- 
range a quantum system at a given excitation energy E. There- 
fore, the entropy of a nuclear system can give information on 
the underlying nuclear structure. The microcanonical entropy 
is given by 



S(E)=k B )nQ.(E), 



(5) 
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1 We chose to apply the old parameterization of | 4] instead of the more recent 
one in | 5] because new (y,n) data | 21] showed that the slope of the radiative 
strength function in 117 Sn Q|] became too steep using the values of Q). 
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FIG. 2: Normalized level density of a) 116 Sn and b) 117 Sn. Results 
from the present work are shown as full squares. The data points 
between the arrows are used for the fitting to known data. The level 
density at lower excitation energy obtained from counting of known 
discrete levels is shown as a solid line, while the level density calcu- 
lated from neutron resonance spacings is shown as an open square. 
The back-shifted Fermi gas level density used for the interpolation is 
displayed as a dashed line. 



where £l(E) is the multiplicity of accessible states and kg is 
Boltzmann's constant, which we will set to unity to give a 
dimensionless entropy. 

The experimental level density p(E) is directly propor- 
tional to the multiplicity £l(E), which can be expressed as 

£l(E)=p(E)-[2(J(E)) + l}. (6) 

Here, (J(E)) is the average spin at excitation energy E and 
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FIG. 3: a) Experimental entropies for 116 ' I17 Sn, and b) the entropy 
difference. By fitting a straight line to the entropy difference in the 
energy region 2.4 — 6.7 MeV, an average value of AS ~ 1.6 kg was 
obtained. 



the factor [2 (J{E)) + 1] thus gives the degeneracy of magnetic 
substates. As the average spin is not well known at all exci- 
tation energies, we choose to omit this factor and define the 
multiplicity as 



Q(E)=p(E)/p , 



(7) 



where the denominator is determined from the fact that the 
ground state of even-even nuclei is a well-ordered system with 
zero entropy. The value of po = 0.135 MeV -1 is obtained 
such that S = lnil ~ for the ground state region of 116 Sn. 
The same po is also applied for 117 Sn. In Fig.[3]the resulting 
entropies of 116 ' 117 Sn are shown. 

The entropy carried by the valence neutron can be estimated 
by assuming that the entropy is an extensive (additive) quan- 
tity l22ll . With this assumption, the experimental single neu- 
tron entropy is given by 



AS = S( 117 Sn)-S( U0 Sn 



116 c 



(8) 



the excitation energy increases, and above 2 — 3 MeV we esti- 
mate the single neutron quasiparticle to carry about AS ~ 1.6 
in units of Boltzmann's constant. This agrees with previous 
findings from the rare-earth region lf22ll . 

Both entropy curves display step-like structures superim- 
posed on the general smooth increasing entropy as a function 
of excitation energy. At these structures, the entropy increases 
abruptly in a small energy interval before it becomes a more 
slowly increasing function. 

The first low-energy bump of 116 Sn is connected to the first 
excited 2 + state at E = 1.29 MeV and the second excited + 
state at E = 1.76 MeV. Similarly, the first bump in the en- 
tropy of 117 Sn is connected to the first excited states in this 
nucleus. The next structures are probable candidates for the 
pair-breaking process. Microscopic calculations based on the 
seniority model indicate that step structures in the level den- 
sity can be explained by the consecutive breaking of nucleon 
Cooper pairs 12311 . 

The bumps present in the Sn level densities are much more 
outstanding than previously measured for other mass regions 
by the Oslo group. One explanation of the clear fingerprints 
could be that since the Z = 50 shell is closed, the break- 
ing of proton pairs are strongly hindered and thus do not 
smooth out the entropy signatures for the neutron pair break- 
ing. Therefore, it is very likely that the structures are due to 
pure neutron-pair break up. 

We have investigated the structures further by introducing 
the microcanonical temperature given by 



dS_ 
dE 



(9) 



By taking the derivative, small changes in the slope of the 
entropy are enhanced. This is easily seen in Fig. |4j where 
the microcanonical temperatures of 116117 Sn are displayed. 
The striking oscillations of the temperature is a thermody- 
namic signature for such small systems as the nucleus. Only a 
few quasiparticles participate in the excitation of the nucleus. 
Since the system is not in contact with a heat bath with a con- 
stant temperature, the system is far from the thermodynamic 
limit. 

In Fig.HJ the bumps below w 2 MeV are connected to the 
low-lying excited states in 116117 Sn. However, the peak-like 
structure centered around E = 2.8 MeV in 116 Sn and around 
E = 2.6 MeV in 1 17 Sn could be a signature of the first break- 
up of a neutron pair. Above E = 4.5 and 3.6 MeV in 116,117 Sn, 
respectively, the temperature appears to be constant on the av- 
erage, indicating a more continuous breaking of further pairs. 



Recently [24], the criticality of low-temperature transitions 
was investigated for rare-earth nuclei. We apply the same 
method here and investigate the probability P of the system 
at fixed temperature T to have excitation energy E, i.e., 

P(E,T)=a(E)esp(-E/T)/Z{T), (10) 

where the canonical partition function is given by 



From Fig. [3] we observe that AS becomes more constant as 



Z(T)= Cl(E')exp(-E'/T)dE f . 



(11) 
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FIG. 4: Microcanonical temperatures of a) 116 Sn, and b) 117 Sn. 



Lee and Kosterlitz showed 1I25I l26ll that the function 
A(E,T) = — ]nP(E,T), for a fixed temperature T in the vicin- 
ity of a critical temperature T c of a structural transition, will 
exhibit a characteristic double-minimum structure at ener- 
gies E\ and Ei. For the critical temperature T c , one finds 
A(E{ , T c ) =A(E2,T C ). It can be easily shown that A is closely 
connected to the Helmholtz free energy, and that this condi- 
tion is equivalent to 



F C (E 1 )=F C {E 2 ), 



(12) 



which can be evaluated directly from our experimental data. 
It should be emphasized that F c is a linearized approximation 
to the Helmholtz free energy at the critical temperature T c ac- 
cording to 



pair. The values found are T c = 0.58(2) and 0.71 (2) MeV for 
116,117 Sn, respectively. Furthermore, we observe a potential 
barrier AF C of about 0.25 MeV between the two minima E\ 
and E2- The potential barrier indicates the free energy needed 
to go from one phase (no pairs broken) to another (one broken 
pair) at the constant, critical temperature T c . 

In the process of breaking additional pairs, the structures 
are expected to be less pronounced. Indeed, in the lowest pan- 
els of Fig. [5] the free energy is rather constant for excitation 
energies above E ss 5.0 and 4.2 MeV for 116 ' 117 Sn, respec- 
tively. Instead of a double-minimum structure, a continuous 
minimum of F c appears for several MeV of excitation ener- 
gies. This demonstrates clearly that the depairing process in 
tin cannot be interpreted as an abrupt structural change typical 
of a first order phase transition. For this process we evaluate 
the critical temperature by a least % 2 fit of Fq to the experimen- 
tal data. The excitation energies (E\ and E2) and temperatures 
for the phase transitions are summarized in Table [TT] 



TABLE II: Phase transition values deduced for 116117 Sn. 



Breaking of 



TTF 

E1-E2 
(MeV) 



Sn 
(MeV) 



Ei-E 2 T c 
(MeV) (MeV) 



F C (E)=E-T C S(E). 



(13) 



one pair 2.4-4.2 0.58(2) 2.3-3.7 0.71(2) 
two or more pairs 5.0-8.0 0.69(8) 4.2-6.0 0.63(4) 



According to the linearized free-energy calculations, the 
first neutron pair is broken for excitation energies between 
2.4-4.2 MeV in 116 Sn. Comparing with the upper panel of 
Fig. |4] we see that this coincides with the excitation-energy 
region where a peak structure is found. This bump is thought 
to represent the first neutron pair break-up as discussed pre- 
viously in the text. The average temperature of the bump is 
0.59(2) MeV, in excellent agreement with the deduced crit- 
ical temperature T c = 0.58(2). Similarly, for 117 Sn in the 
lower panel of Fig. |4j we find that a peak-like structure with 
an average temperature of 0.74(2) MeV is present between 
E = 2.3 — 3.7 MeV. Compared to the critical temperature 
T c = 0.71(2) MeV, the values agree satisfactory. 

For the continuous break-up process characterized by a zero 
potential barrier (AF C « MeV), we estimate for 116 Sn an 
average microcanonical temperature of 0.75(6) MeV in the 
excitation-energy region E = 5.0 — 8.0 MeV, in reasonable 
agreement with T c = 0.69(8) calculated from the linearized 
Helmholtz free energy. In the case of 117 Sn, the average 
microcanonical temperature for excitation energies between 
4.2 - 6.0 MeV is found to be 0.64(3) MeV, which agrees very 
well with T c = 0.63(4) MeV. This gives further confidence in 
our interpretation of the data as two independent methods give 
very similar results. 



Linearized free energies F c for certain temperatures T c are 
displayed in Fig. In the upper panels, 116117 Sn data are 
shown where the condition F C (E\) — F C (E2) = F$ is fulfilled. 
Each nucleus shows a double-minimum structure, which we 
interpret as the critical temperatures for breaking one neutron 



V. CONCLUSIONS 

The nuclear level densities for 116117 Sn are extracted from 
particle-y coincidence measurements. New experimental re- 
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FIG. 5: Linearized Helmholtz free energy F c at critical temperature T c for a) 116 Sn, and b) 117 Sn displaying a characteristic double-minimum 
structure. The constant level Fq, which is connecting the minima, is indicated by horizontal lines. A continuous minimum of F c is shown for 
c) 116 Sn,andd) 117 Sn. 



suits are reported for the level density in 116 < 117 Sn for excita- 
tion energies above 1.5 MeV and 3.5 MeV up to S„ — 1 MeV. 
The level densities for both nuclei display prominent step- 
like structures. The structures have been further investigated 
by means of thermodynamic considerations: microcanonical 
entropy and temperature, and through calculations of the lin- 
earized Helmholtz free energy. Both methods give consistent 
results, in strong favor of the pair-breaking process as an ex- 
planation of the structures. 
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